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ABSTRACT 

In the standard picture of disc galaxy formation, baryons and dark matter receive the 
same tidal torques, and therefore approximately the same initial specific angular mo- 
mentum. However, observations indicate that disc galaxies typically have only about 
half as much specific angular momentum as their dark matter haloes. We argue this 
does not necessarily imply that baryons lose this much specific angular momentum as 
they form galaxies. It may instead indicate that galaxies are most directly related to 
the inner regions of their host haloes, as may be expected in a scenario where baryons 
in the inner parts of haloes collapse first. A limiting case is examined under the ide- 
alised assumption of perfect angular momentum conservation. Namely, we determine 
the density contrast A, with respect to the critical density of the Universe, by which 
dark matter haloes need to be defined in order to have the same average specific 
angular momentum as the galaxies they host. Under the assumption that galaxies 
are related to haloes via their characteristic rotation velocities, the necessary A is 
~ 600. This A corresponds to an average halo radius and mass which are ~ 60% and 
~ 75%, respectively, of the virial values (i.e., for A = 200). We refer to this radius as 
the radius of baryonic collapse Rbc, since if specific angular momentum is conserved 
perfectly, baryons would come from within it. It is not likely a simple step function 
due to the complex gastrophysics involved, therefore we regard it as an effective ra- 
dius. In summary, the difference between the predicted initial and the observed final 
specific angular momentum of galaxies, which is conventionally attributed solely to 
angular momentum loss, can more naturally be explained by a preference for collapse 
of baryons within Rbc, with possibly some later angular momentum transfer. 

Key words: galaxies - formation, galaxies - evolution, galaxies - kinematics and 
dynamics, galaxies - fundamental properties. 
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C3 1 INTRODUCTION 



In the standard picture of disc galaxy formation (e.g., Fall| 
fc Efstathiou||1980| |Dalcanton, Spergel, fc Summers|[l997| 



Mo, Mao, & White 19981, galaxies consist of a dissipa- 



tive baryonic component and a non-dissipative dark matter 
component. Galaxies form hierarchially, and in this process, 
baryons and dark matter acquire the same specific angu- 
lar momentum (j) via tidal-torques. This is because tidal- 
torques are most effective in the linear and the trans-linear 
regimes, when baryons and dark matter are well-mixed. 
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The dark matter then collapses non-dissipatively, and the 
baryons dissipatively, likely with some cloud-cloud collisions 
and possibly shocks (processes which are expected to re- 
arrange j but not remove it). The baryons form rotating 
centrifugally-supported discs at the centres of the potential 
wells. For a review of this scenario see Fall (2002 1. This 
standard picture is able to correctly predict galaxy prop- 
erties such as scale-lengths and sizes if the baryons retain 
most of their initial j. It has been extended to include ad- 
ditional physics effects and larger samples of galaxies by 
e.g., I White & Frenkl (l99ll), ICole et all Jl994|, ISomerville 



& Primack (19991, de Jong & Lacey (2000), Van den Bosch 



(2001), Hatton et al. (20031, and Dutton (2009) 



In order for this scenario to correctly predict galaxy 
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properties, the baryons must retain a large fraction of their 
initial angular momentum. However, early numerical sim- 
ulations of galaxy formation contradicted this expectation 



(Katz&Gunnl991 Navarro & Benz 1991 Navarro & White 



19941. They found a factor of ~ 30 loss of angular momen- 



tum for simulated galaxies, and referred to this as an "angu- 
lar momentum catastrophe." As simulations improved over 
the years, it became clear that much of this catastrophe was 
actually a numerical artifact: too little resolution and too 
much numerical viscosity (see e.g., Governato et al. 2010 
Brooks et aLl[20TT| |Keres et aLl|2011| |Brook et al. fboil 



Kimm et al. 



2011 



and references therein). Another possi- 
ble contribution to solving the angular momentum problem 
may be through feedback effects which can delay baryons 



from falling onto discs (e.g., Weil, Eke, fc Efstathiou||1998[ 
Sommer-Larscn, Gclato, & Vedel 1999; Eke. Efstathiou, & 



Wright 2000 Thacker & Couchman 20011. With high nu- 



merical resolution and some feedback, galaxy simulations 
are now at a stage where angular momentum loss may be a 
relatively minor problem. In this paper, we explore another 
option: that the discs of galaxies draw baryons mainly from 
the inner parts of dark matter haloes. Some of the baryons 
in the outer parts may have not yet collapsed onto the discs. 

The angular momentum catastrophe prompted compar- 
isons of the j of simulated haloes to that of observed galax- 
ies. In these studies, the j of dark matter haloes is measured 
out to the virial radius, Rvir, which is standardly defined 
as -Ra=200, and is the effective radius at which the dark 
matter ceases to collapse into the ha lo. |Navarro fc Stein-| 



mctz 



(20001 and Burkert & D'Onghia (2004) found that ob- 



served galaxies have 45% and 70% of the j of their expected 
host haloes in simulations, respectively, under the assump- 
tions that galaxies can be related to simulated host haloes 
via characteristic rotation velocities directly and via a scal- 
ing factor, respectively. Recently, |Dutton fc van den Bosch] 
(20121 found that the spin parameters of observed galaxies 
are ~ 60% of those of simulated haloes. These studies are 
consistent once differences in assumptions and approxima- 
tions are accounted for. 

Studies which compare the total j predicted for haloes 
by numerical simulations to that observed for galaxies all 
assume that the effective outer halo radius from which the 
baryons collapse (defined here as Rbc) is equal to Rvir- Be- 
cause baryons in the inner parts of haloes will have higher 
cooling rates and more frequent cloud-cloud collisions, it is 
reasonable to expect that they form the galaxies, and that 
baryons from larger radii are not captured. Although Rvir 
has traditionally been identified with Rbc, these two radii 
are governed by different physics (dissipative versus non- 



dissipative), and need not be related, as emphasized by Fall 



(2002|. The only requirement is that Rbc must be interior 
to Rvir, since baryons cannot collapse from unvirialized re- 
gions. The purpose of this paper is to determine the effect 
of relaxing the assumption that Rvir and Rbc are equal on 
the difference in j between galaxies and haloes. We assume 
for simplicity that the boundary between the collapsed and 
uncollapsed baryons is a sharp one. In reality, it will be a 
gradual boundary because some of the baryons in the halo 
within Rbc might not collapse, and some baryons outside 
of Rbc might. Therefore, we regard Rbc as the effective 
boundary between these two regions. 

In this paper, we ask the following question: If galaxies 



formed from all the baryons in haloes out to Rbc, and be- 
yond this radius the baryons remained in the halo, what is 
the value of Rbc required to match the j of galaxies? We 
address this question by comparing the j observed for disc 
galaxies with that of their expected dark matter haloes mea- 
sured within a range of halo radii. For disc galaxies, j can 
be measured from observations of surface brightness profiles 
and rotation curves. For dark matter haloes, we must resort 
to numerical simulations. 

This paper is organised as follows: In §2, we measure 
j of dark matter haloes in a cosmological dark matter-only 
simulation. We investigate its dependence on the halo ra- 
dius within which j is measured and the halo radius at 
which the rotation velocity is measured. The resulting pre- 
dictions of dark matter halo j are compared to j measured 
for a large observational sample of local galaxies for which 
the completeness is known in §3. A discussion of the re- 
sults is in §4. We adopt a ACDM concordance universe 
(On = 0.24, fl A = 0.76, ft = H / [100 km s" 1 Mpc -1 ] = 
0.73, as = 0.77, n — 0.958), i.e., within one standard devia- 



tion of both the WMAP 3 and 5 year best estimates ( Spergel 
|et al.|2007||bunkley et al.|2009| ). All logarithms are base ten. 



2 N-BODY SIMULATION OF DARK MATTER 
HALOES 

To quantify the dependence of dark matter halo j on how 
the outer radius of a halo is defined, we look to a suite 
of cosmological N-body simulations of dark matter haloes. 
These simulations include only dark matter and gravity (i.e., 
neither baryons nor hydrodynamics). As discussed in the 
Introduction, if the baryons in a given dark matter halo are 
initially distributed in the same manner as the dark matter, 
and they later cool to form a disc while conserving j, then 
the j of the galaxy should be equal to that of the virialized 
region of the dark matter halo. However, if baryons collapse 
progressively from the inner to the outer parts of haloes, and 
they have not finished collapsing (or, if some baryons never 
collapse), then galaxy j may be expected to reflect that of 
dark matter haloes within a given radius, Rbc- 

To predict the distribution of j among dark matter 
haloes, a large N-body simulation is needed which can model 
the acquisition of angular momentum for even the slow- 
est rotating galaxies in our sample (125 km s" 1 ; §3). The 
simulation we adopt is part of the Horizon Project suite 
(http://www.projet-horizon.fr). This follows the evolution 
of a cubic cosmological volume of 100 h" 1 Mpc on a side 
(comoving) containing ~ 134 million dark matter particles 
(512 3 ). It starts at z — 99 and is evolved using the publicly 
available treecode Gadget 2 ( Springcl 2005 ) with a softening 
length of 5 h _1 kpc (co-moving). The adopted cosmology re- 
sults in a dark matter particle mass of 6.83 x 10 s M©. Dark 
matter haloes and the subhaloes they contain are identified 



with the AdaptaHOP algorithm ( Aubert et al.|2004). The halo 



centres are positioned on the densest dark matter particle 



located in the most massive substructure (see Tweed et al. 
2009 for details). The total number of haloes and subhaloes 
in the simulation volume at z — with more than 100 par- 
ticles within -R200 and with circular velocities at this radius 
which are greater than 100 km s _1 is 9661. 

The j of a halo is measured within a range of radii 
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Figure 1. For simulated dark matter haloes at z 0, the relations between ja (for A = 200, 2000, and 20,000) and rotation velocities 
V200 and V^o.OOO are shown. Individual haloes are plotted as grey points, binned averages are shown as black triangles, and the rms 
scatter is shown as black error bars. Contours in volume density are shown for 2 and 20 X 10 — 5 haloes per 0.1 in log j'a and per 0.1 log 
V, per Mpc 3 . The shapes of the distributions are similar for j'a whether V200 or ^20,000 is adopted. As A increases, the normalisation 
of the relation between ja and V decreases, but the slope and scatter do not change greatly. Similar relations are found for V2000, but 
are not shown to avoid redundancy. 



as follows. First, the halo is divided into 100 radial ellip- 
soidal shells, where the axis ratios of the ellipsoid are ob- 
tained by computing the inertial tensor of all the parti- 
cles in the halo. Halo circular radii are defined as the cube 
root of the radii of the three major axes of each ellipsoid. 
Next, the vector angular momentum of the particles in each 
shell is calculated, and the angular momenta of the shells 
is summed vectorially from the inner-most shell to the radii 
specified before taking its modulus. The mass of a halo is 
measured in an analogous manner, and j is simply the an- 
gular momentum divided by the mass within a given ra- 
dius. Selected radii, Ra, are defined by the density of the 
haloes with respect to the critical density of the universe 
(A = p(r < i?A)/pcrit)- Specific angular momenta measured 
within these radii are defined as ja- Circular velocities at 
these radii are Va = (GM a /Ra) 1/2 , where M A and Ra 
are the mass and radius of the halo defined by A, and G is 
the gravitational constant. The ranges of A, RA/R200, and 
Ma/M 20 o probed are 50-20,000, 1.70-0.09, and 1.24-0.13, 
respectively. 



In Figure [T] relations between halo ja and Va are 



shown|jHalo j is measured within -R200, ^2000, and i?20,ooo, 
and halo V is measured at -R200 and ifeo.ooo- We do not show 
results for halo V measured at -R2000 since they do not differ 
significantly from those for i?200 or 7?2o,ooo- The radii -R2000 
and i?20,ooo correspond to 34% and 9% of -R200, respectively, 
on average Only haloes with more than 100 particles are 
retained, except for measurements of ^20,000 and V2o,ooo for 
which haloes with more than 50 particles are used. For these 
50-particle haloes, the intrinsic relations remain the same, 
but the scatter is increased slightly due to increased Poisson 
noise. The shapes of all the distributions are similar in terms 
of slope and scatter, and are therefore approximately inde- 
pendent of the radius for which j or V are measured. The 
slope flattens slightly with increasing A, and the scatter re- 
mains about the same. We will quantify this in the following 



1 There is a drawback to a plot of j versus V, namely both axes 
incorporate factors of V, and a relation is expected by construc- 
tion (e.g., Freeman 1970). Because the local relation between 
galaxy V and stellar mass is tight (e.g., |Bell fc de Jong|[2001| 
|Kassin, de Jong, fc Weiner|20u6] ), there is a similarly tight rela- 
tion between between j and stellar mass (e.g., |Fall|1983l l, which 
is not expected by construction. 
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Figure 2. These plots are the same as in Figure^for V20O1 except here observed galaxies are also shown (Vfi a t of the galaxies is adopted 
as a characteristic rotation velocity and is plotted on the horizontal axis). Individual galaxies are plotted as green points, and contours 
in volume density are shown in blue for 2 and 20 X 10 — 5 galaxies per 0.1 in log j/± and per 0.1 in log V200 P er Mpc 3 . Binned averages 
for the galaxies are shown as blue triangles, and the rms scatter is denoted by error bars. The scatter in j for galaxies is about half of 
that of the haloes. Under the assumption that characteristic galaxy and halo rotation velocities are equal (i.e., Vfi a t = ^200)1 galaxies 
have on average a factor of ~ 2 less j than haloes defined with A = 200, a factor of ~ 2 more j than haloes defined with A = 2000, and 
a factor of ~ 5.6 more j than haloes defined with A = 20, 000. 



section. However, the normalisation is strongly dependent 
on A: it decreases by factors of ~ 3 and ~ 6 for 10 and 
100-fold increases in A, respectively. A decreasing normal- 
isation with increasing A is a consequence of how angular 
momentum is distributed in galactic haloes, with most of 
the angular momentum located in the outer parts. As we 
increase A, we exclude more and more of the outer parts of 
the haloes, and the angular momenta decrease, as illustrated 
by the simple analytic treatment in |Fall| ( |1983| Section 4) . 
In this paper, we quantify this decrease more precisely using 
numerical simulations. 



3 COMPARISON WITH OBSERVATIONS OF 
DISC GALAXIES 

The goal of this section is to place measurements for disc 
galaxies on Figure[l] To do so, we need to (1) adopt a galaxy 
sample for which the completeness is well-defined and which 
has the necessary data available to derive circular velocities 
and j, and (2) relate galaxies to simulated host dark matter 
haloes. 

To address the first need, a large sample of 456 galax- 



ies from Mathewson, Ford, & Buchhorn ( 1992 1 and com- 
pleteness ineasuTelnlmtsHfrOTn^de^ ( |2000| ) are 
adopted. Details of this sample are given below. The large 
size of and the data available for the sample necessitates 
simple estimates of j. Therefore, we estimate j as 2Vn a ,tr<i, 
where Vfl a t is the rotation velocity on the flat part of the ro- 
tation curve and Td is the scale-length of the galaxy disc. 
This approximation is exact for an exponential disc and 
a flat rotation curve. Uncertainties in estimates of j are 
~ 15%, which are dominated by errors in measurements 
of r s (mainly due to errors in sky background subtraction) 
and galaxy distances. There are two minor effects on esti- 
mates of j, which we do not take into account, but which 
work in opposite directions. On the one hand, galaxies have 
rising rotation curves in their centres, and this causes the 



formula to slightly overestimate j. On the other hand, most 
galaxies are expected to have extended gas discs, but with 
very little mass, which would cause the formula to slightly 
underestimate j. 

The galaxy sample used is a sub-sample of the ESO- 
Uppsala Catalog of Galaxies (Lauberts: 1982) which was se- 
lected by eye from photographic plates. It is only incomplete 
for very late Hubble types (T > 6, i.e., later than Scd; de 
Jong fc Laceyl|2000| . Values of V fl at were determined from 
optical and radio observations. For the optical data, V% at 
was defined as half the difference between the maximum and 
minimum velocities of the Ha rotation curves. For the radio 
data, Vkat was defined as half the width of the HI profile be- 
tween points where the intensity falls to 50% of the highest 
values; these values were then corrected for dispersion and 
converted to optical rotation velocities by multiplying by 
1.03 and then subtracting 11 km s _1 (see §3.4 and Figure 



5 of |Mathewson, Ford, fe Buchhorn||1992| ). Disc half light 
radii, which are the result of J-band bulge-disc decomposi- 



tions from de Jong & Lacey (2000), are converted to disc 



scale lengths by dividing by 1.679 (the exact ratio of the 
half-mass radius to the scale radius for a pure exponential 
disc). Only those galaxies with rotation velocities greater 
than 125 km s _1 are used. This helps us to avoid galaxies 
with rotation curves which do not flatten out at the radii 
measured. The distribution of galaxies in j versus Va a t does 
not differ significantly from the galaxy sample commonly 



used in the literature (Courteau et al. 20071, but it has a 
better completeness. 

To address the second need, and relate galaxies to the 
dark matter haloes in Figure [l] we assume for simplicity that 
the characteristic rotation velocity of a galaxy (which we 
take to be Vfiat) and that of its host halo at R200 are equal. 
For a massless disc in a |Navarro, FrenkTfc White ( 1996 ) halo, 
V c at the location of the galaxy can be about half its value 
at i?200- However, the self-gravity of the baryons is expected 
to increase V c in the inner parts of haloes. The amount by 
which it increases is difficult to calculate theoretically, so we 
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Figure 3. The average difference between galaxy and halo log j, < log j ga i ax i cs - log j'a >, is shown as a function of A, -Ra/^200! an d 
M a /M 200 , in panels a, b, and c, respectively. Points demarcate discreet values, and solid lines simply connect the points. There is no offset 
between galaxy and halo log j for A = 578^,, which corresponds to Rbc = Ra=57s/^200 = 0.63^q'q^, and M&—§j%/M2oo = 0.74±0.1. 



look to observations. |Dutton et al.| ( |2010[ ) find a very small 
conversion factor between V c at -R200 and at the location 



of galaxy discs. In their analysis, Dutton et al. (2010) com- 



bined dark halo masses measured from satellite kinematics 
and weak gravitational lensing to show that V2.2 — V200 for 
V2.2 = 90 — 260 km s _1 , where V2.2 is the galaxy rotation ve- 
locity measured at 2.2 /-band scale lengths. This equivalence 
is also consistent with semi-analytic models of galaxy forma- 
tion which require a similar ratio between galaxy and halo 
velocities to simultaneously match the local Tully-Fisher re- 



lation and galaxy luminosity function (e.g., Dutton & van 
|den B osch 2009] and references therein). 

In Figure |2j we compare the distribution of j versus 
Vk a t for galaxies described in this section with the distribu- 
tions of ji'200, .72000, and .720,000 versus V200 for dark matter 
haloes from Figure [T] As discussed above, it is assumed that 
haloes have the same rotation velocities as the galaxies they 
host, so they can be directly compared in Figure [2] The 
halo relations from Figure [l] for V-jo.ooo are not shown be- 
cause they are not significantly different from those for Vioo- 
We fit a linear relation to the galaxies using 100 bootstrap 
re-samplings and a generalised least squares fitting routine 
I Weiner et aT|2006b[ ), which gives a slope of 2.5 ± 0.1 rms. 
We also fit a linear relation to the haloes in Figure [2] for j'200 
versus V200 for circular velocities which span the velocity 
range of the galaxies, 125 < V200 < 315 km s . This results 
in a slope of 1.92 ±0.02 rms. The distribution of galaxies has 
a similar slope to that of the haloes, as fo und by|Fall|(|1983[) 
and others (e.g., |Mo, Mao, fc White|1998||Navarro fc Stein- 
metz|2 000), and approximately half the average rms scatter 
(0.15 dex versus 0.27 dex). The lower scatter compared to 
the haloes is related to the finding by de Jong & Lacey ( 2000 1 



that the width of the observed scale-radius distribution of 
galactic discs is narrower than that expected from the distri- 
butions of halo spin parameters in cosmological simulations. 
For the halo .72000 versus V200 and ,720,000 versus V200 rela- 
tions, the slopes are 1.80T0.02 and 1.26T0.02, respectively, 
and the average rms scatters are 0.28 and 0.26, respectively. 
The slopes flatten slightly with increasing A, but the scatter 
remains constant to within errors. Given all the factors not 
included in our simple picture, we consider it remarkable 



how similar the galaxy and halo slopes are. The main result 
of this paper is encapsulated in the much larger difference 
in normalisation between galaxies and haloes. We choose to 
measure this difference at approximately the center of the 
distributions, at log K t=2.35 (V ro t = 224 km s" 1 ). The 
average normalisation of the galaxies is less than that of 
the haloes for j'200 by a factor of ~ 2 (0.30 dex), consistent 
with previous studies (e.g., Navarro fc Steinmetz|2000" Dut- 



ton & van den Bosch 20121. The average normalisation of 



the galaxies is greater than that of the haloes for j'2000 and 
.720,000 by factors of ~ 2 (0.30 dex) and ~ 5.6 (0.75 dex), 
respectively. 



We quantify the dependence of j'a on A as follows. We 
start by measuring halo j for a range of A and compare them 
with j measured for the galaxy sample, as in Figure [2] In 
Figure [3] we show the average difference between galaxy and 
halo j (measured at log V ro t = 2.35) as a function of A, halo 
outer radius in terms of -R200 (i-e., Ra/ '-R200), and halo mass 
in terms of M 2 oo (i-e., Ma/M 2 oo)- The quantities 7?a/ '-R200 
and A/a /M200 are average values for all haloes with circular 
velocities which span 125 < V200 < 315 km s" 1 . The value 
of A at which the average j of galaxies and haloes match 
(i.e., < log jg a i a xi cs - log j A > = 0) is 578^31- The halo j 5 78 
versus V 2 oo relation has a slope of 1.89T0.03 and an average 
rms scatter of 0.28 dex over the velocity range of the galax- 
ies. This value of A corresponds to Rbc = Z?a=57s/^200 = 
0.63lo.i and Ma=578/M 20 o = 0.74 ±0.1. We calculate these 
values by interpolating the curves in Figure |3j and the er- 
rors by considering the limiting case that galaxy j values 
are systematically over and under-estimated by the assumed 
measurement uncertainty. If baryons conserved j perfectly 
during galaxy formation, then the collapse radius Rbc is 
~ 63% of the virial radius i?200- This portion of the haloes 
contains on average 74% of their mass, and if baryons and 
dark matter are initially well-mixed, the same percentage 
of the baryons. However, as discussed in the next section, 
this radius and mass fraction are probably not simple step 
functions; therefore we regard them as "effective" quantities. 
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4 DISCUSSION 

In this paper, we determine the extent to which the ap- 
proximate factor of 2 discrepancy between the j of galaxies 
and their expected host dark matter haloes is sensitive to 
the conventional assumption that Rbc = i?200- This differ- 
ence in j is usually attributed to loss of baryonic j during 
galaxy formation. However, there is no physical reason for 
the assumption that these radii are equal to at least within 
a factor of ~ 2, as emphasized by Fall (20021. This is be- 



cause different physics governs each, namely dissipational 
and dissipationless physics for Rbc and Rvir, respectively. 
The only constraint on the relationship between these radii 
is that Rbc must be interior to Rvir since baryons cannot 
collapse from a region that is not incorporated into the halo. 
A Rbc which is interior to Rvir is a natural expectation in 
the standard theory of galaxy formation where the inner 
parts of haloes collapse first. As Rbc decreases, the discrep- 
ancy between the j of galaxies and haloes is alleviated. We 
show that the discrepancy can be explained entirely by a 
Rbc which is ~ 60% of Rvir- 

To do so, we determine the value of Rbc at which the j 
of galaxies and haloes match. This is done by comparing the 
distribution of j observed for a sample of local disc galaxies, 
for which the completeness is understood, to that predicted 
for their host dark matter haloes from a dark matter-only 
simulation of the Universe. It is assumed that galaxies and 
haloes can be related directly via their rotation velocities. 
The necessary value of the density contrast A needed to 
define the haloes which have the same average j as galaxies 
is ~ 600. This corresponds to an average effective Rbc which 
is ~ 60% of -R200, and an average halo mass which is ~ 75% 
of A/200- Therefore, if galaxies formed from baryons initially 
present in the inner parts of their host haloes and conserved 
j perfectly, the baryons would come from within Rbc and 
would comprise this percentage of the baryons in the halo. 

Even under the assumption of perfect conservation of j, 
Rbc is not likely a sharp boundary. The baryons which form 
the galaxy may only on average come from within Rbc, with 
most material originating from smaller radii, but some from 
more distant radii. In addition, the smaller scatter of the 
galaxies in j versus V compared to that of the haloes may 
indicate a mechanism by which only selected baryons form 
the disc, regulatory processes which act upon the baryons, 
and/or haloes which form non-disc galaxies. This is because, 
in our simple picture, the initial distribution of baryons in 
j versus V is expected to mirror that of the dark matter. 
Therefore, if only selected baryons formed discs or regula- 
tory processes acted upon them during disc formation, it 
may be expected that the baryons which form the discs 
would have a narrower distribution in j versus V. In ad- 
dition, since we compare the predicted properties of dark 
matter haloes with those of disc galaxies, not ellipicals which 
rotate slower than discs, it stands to reason that the com- 
bined population of discs and ellipticals would be broader 
in j versus V ( |Fall|1983| ). 

Eventually, it should be possible to compute Rbc from 
hydrodynamical and dark matter simulations of galaxy for- 
mation in a cosmological context. Current simulations may 
have spatial and mass resolutions that are too coarse to 
model accurately the complex processes expected to be at 
play, such as gas shocks, cloud-cloud collisions, and a mul- 



tiphase medium. These processes affect the rate at which 
the baryons collapse, but they have may relatively little in- 
fluence on the angular momentum of the resulting galactic 
discs. 

A number of phenomena can alter the j of galaxies (see 
|Fall||2002| and |Romanowsky fc Fall||2012| for more complete 
discussions of these phenomena). For example, torques ex- 
erted between the dark matter and the baryons could in 
principle spin up the halo and spin down the disc. Minor 
mergers might also affect the j of galaxies. In addition, feed- 
back from star-formation can alter j differently depending 
on how it varies with radius. Material in outflows may be 
launched from inner or outer radii, or both. If material is 
primarily removed from the inner or outer parts of galaxies, 
galaxy j will increase or decrease, respectively. If feedback 
is active but independent of radius, then there would be no 
change in j. We expect some of these phenomena to alter 
the j of discs, but whether they have a major or a minor 
effect on galaxy j is still uncertain. In order to perform a 
more detailed comparison of galaxies and haloes, we need a 
better understanding of the processes of j transfer in galaxy 
formation, and whether outflows can change the j of galax- 
ies. 

In summary, the difference between the predicted initial 
and the observed final j of galaxies, which is conventionally 
attributed solely to angular momentum loss, hinges on the 
loosely motivated assumption that all the baryons within 
Rvir collapse to form galaxies. There is no physical reason 
why this has to be the case. If baryons in the inner parts of 
haloes collapse first, as is expected, then the j discrepancy 
between galaxies and haloes can be fully explained by a col- 
lapse radius Rbc which is ~ 60% of the virial radius Rvir- 
In the future, baryons from progressively larger radii in the 
halo may collapse, and at some point in time Rbc might 
equal Rvir- In reality, it may be that a combination of a 
preference of collapse of the inner parts and some j trans- 
fer between baryons and dark matter is needed to solve the 
problem. 
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